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ABSTRACT 

Two types of integrable coupled nonlinear Schrodinger (NLS) equations are derived by 
using Zakharov-Shabat (ZS) dressing method. The Lax pairs for the coupled NLS 
equations are also investigated using the ZS dressing method. These give new types of the 
integrable coupled NLS equations with certain additional terms. Then, the exact solutions 
of the new types are obtained. We find that the solution of these new types do not always 
produce a soliton solution even they are the kind of the integrable NLS equations. 



1 Introduction 

One of the most remarkable discoveries in soliton the- 
ory is the inverse scattering method (ISM). It is well- 
known that there are several types of the coupled 
nonlinear Schrodinger (NLS) equations whose com- an( j 
plete integrability can prove by some methods in the 
ISM such as AKNS method §J, Wadati method, 
ZS dressing method, and so on [2 — 5]. In this paper, ^ ^ ^ 

we introduce two new types of the integrable coupled e secon type, 
NLS equations derived using ZS dressing method. 
The first type of the equations is 
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iu t + xu 



iw t 



(1.1) 



(1.2b) 
(1.2c) 

(1.2d) 

Qi 

Q2, 

(1.3) 



is a coupled NLS equation with complicated addi- 
where the perturbative terms i?i and R 2 , and the tional terms (Qi and Q 2 ) defined as follows 
real parameters and [i are, respectively defined as 
follows 

Rx = -^-s-, (1.2a) 



14; 



and 



* e-mail: elimhi@hotmail.com 



Ql2 + Ql3 
Q2 = Q21 + Q22 + Q23 + Q24 



(1.4a) 
(1.4b) 
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where terms appeared in the related equation provided by 

. . Christodoulides et al. into the complicated terms (Qi 

[\u\ 2 + \w\ 2 j t ^ andg 2 ). 

(1.5a) The present paper consists of the following. In sec- 
tion 2, we perform the ZS dressing method for the 
general NLS equations. In section 3, we propose and 
derive two types of the integrable coupled NLS equa- 
(1.5b) tions using the ZS method. In section 4, we solve the 
solutions of the coupled NLS equations. The last sec- 
tion, section 5, is devoted to the concluding remarks. 

(1.5c) 

2 ZS Dressing Method for the 
General Coupled NLS Equa- 
tions 

2.1 The Integral Operators of The ZS 
Dressing Method 

(1.5e) In 1974, Zakharov and Shabat generalised the Lax 
method using their ZS dressing method. In this sec- 
tion, we will review how the ZS dressing method 
works on the general NLS equations and will also 
(1.5f) show the connection with the work of Lax soon, but 
first we introduce three integral operators [9]. Let, 
in general, F(x,z) and k±(x,z) be NxN matrices 
where 

' ((|«| 2 + H 2 )) 2 \ k+(x,z)=0, H z<x, (2.1a) 

1 " k_ (x,z) = 0, if z > x, (2.1b) 




(1.5d) 



Q 



3s 2 X 



V 



and let ib (x) be an TV- vector. The integral operators 
<&f and $± on ib are defined by 



(i + s (m 2 + m 2 )) 

(l-5g) 

The parameter s is an effective saturation parameter. 
If we put s — > 0, then the system (1.3) reduces to the °9 

integrable Manakov equation [6] : $ F (ip) = / F (x , z) ip (z) dz (2.2) 

iu t + x u ^ T 2/U ( |u| + |iu| 1 u = 

iw t + X w xx T (l u | 2 + \ w \ 2 ^j w = 0, (1.6) 

(ib) — / k+ (x z) ib (z) dz (2 3) 

Equation (1.3) can be identified as the similar ± v ^ ' J ± v > J ^ \ ) > \ ■ ) 

equation with that which has been derived by 

Christodoulides et al. [7], and it is described by a gQ ^ ^ 

system of two coupled nonlinear equations for the °° 

normalized beam envelopes, u(x,t) and w(x,t) [8]. $ + = / k + (x, z)dz, (2.4a) 

However, in eq.(1.3), we have extended the additional J 



for all integrable ip, similarly 
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and = k+ (x, z)+ F (x, z) 

/OO 
k-(x,z)dz. (2.4b) f , N , x 

+ / k+(x,y)F(y,z)dy, for z > x . (2.9) 

-oo J 

t 

We suppose that $_f and <I>± are related by the op- 
erator identity Equation (2.9) is the matrix Marchenko equation for 

k + (x, z). Similarly, if we consider z < x in eq.(2.9), 

(I + <f>+) (7 + $f) = (I + , (2.5) it can be shown that 

oc 

where we assume that (I + <I> + ) is invertible so that f 

k- (x, z)= F (x, z)+ k+ (x, y) F (y, z) dy , (2.10) 

(7 + & F ) = (7 + $+) _1 (7 + $_), (2.6) t 

i.e. the operator (7 + $ F ) is factorisable and 7 ,as which defines fc_ (x, z) in terms of k+ and F. At this 

usual, is the unit matrix. stage we have not restricted the choice of F, and this 

The identity (2.5), on ip, can be written as we will next do for some NLS equations. 

OO OO 

f 7 , x , / x 7 f r,/ \ , / x 7 2.2 The differential Operators of the 

7 k+ (X ' Z) ^^ dz+ J F (»> * W d2 ZS Dressing Method 

£ — OO 

oc / oc \ In common with our previous analyses, we extend the 

+ / k + (x, z) I / 7 1 [z, y) ip (y) dy ] dz definitions of the matrices k± and F, and the vector 

J \ J I ip, so that they all may now depend upon auxiliary 

X \— OO / 

x variable e.g. t, y. We shall describe the evolution of 

/, . . , . . , .„ „. k± and F (in t,y), and hence relate them to certain 

k-(x,z)ip(z)dz, (2.7 , .. y I.'' , . , , . . 

v ' evolution equations by introducing appropriate (lin- 

1 ear) differential oprators. We define the NxN matrix 

and if we require that ip (z) = for z < x then the differential operator Aq \ i = 1, 2 on ip (x; t, y) which 

right-hand side is zero. Furthermore, the double in- has only constant coefficients and which commutes 

tegral may be expressed as with the integral operator <J> F , i-e. 



OO oc 



a£ j $ f - $ F A< !) = 0. (2.11) 

k+(x,y)F(y,z)ip(z)dydz, 

Note that in term $pA^ , operates on ip (x; t, y) 
where the 'dummy' variables have been relabelled by first and thcn this is evaluated on x = z for the appli- 
interchanging y and z. Thus equation (2.7) becomes cation of the operator & F . Further we introduce an 

associated differential operator, A w , which is defined 
by the operator identity 



-CO X 



OO 

= J [k+ (x, z) + F (x, z)] ip (z) dz 



A« (I + $+) = (I + $+) A , 1 = 1,2. (2.12) 

°° It can be shown that eq.(2.12) also holds if $+ is 

/ (k+(x,y)F(y,z))iP(z)dy, (2.8) replaced by $_: 

Theorem 1. If operator (7 + $ F ) commutes with 

, ; . . ,_ , . . ., , , „ „, r the differential operator A^\ and is invertible then 

(note the use here of eq.(2.1a) for all ip (z)). lherefore u (i) 

both its Volterra factors (7 + $±) transform Aq 

into one and the same operator AW . 
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Proof (of theoreml) . The associated operator AW is now obtained from 

eq.(2.12) (for i = 1), 

A« = (/ + * + ) A« (/ + 

((/ + *_)(/ + 1 A^ I V (x; t) + y k+ (x, z- 1) V- (z; t) dz 



= I (iaipt - tb xx ) 

OG 

/" f d d 2 \ 

+ k + (x, z; t) I ia— - ) V (2; *) dz. 



:(A« (/ + $ f )(7 + $_)- 1 ) 
- (7 + $_) A$°(J + (2.13) 

The operator is sometimes referred to as 

'undressed', and AM as the 'dressed' operator. 

Before the general development, we shall consider (2-19) 
an example which will illuminate the meaning of Agaillj we intcgrate by parts to find 
equations (2.11) and (2.12). Let 

A (D = j (ia<L _ jfL) (2 14) / k +^** dz = -k+^ + k+rf + J k+.rfdz, (2.20) 

\ dt dx 2 J ' ' x x 

where a is an arbitrary real value, and I is the where k + = k+(x,x;t), and we assume that fc+, 

3x3 unit matrix. Thus cq.(2.11), when operated on k +z — ► as z — > +00. It is now convenient to set 
tb (x; t), becomes ^ 

^ V h AW=AW+^, 2 = 1,2 , (2.21) 

OO 

= fia°. - ^) [ F (x, z; t) * (z; t) dz S ° that eq '( 2 - 19 ) beCOmCS ( for inStanCC ' 1 = ^ 



J F{x,z;t) (iajt - ip(z;t)dz. 



\'t ! (• + J k + ibdz J + ia J [k +t ib { k i',)d: 



(2.15) 

After integration by parts, it follows that °j? ^ 



00 00 

J F*b zz dz= J F zz ^dz, (2.16) 



OC- 
X 

OO 

= ia J k + tb t dz + k + tb x ~k +z ip~ip — k + , (2.22) 



or 



-co 



for all bounded continuously twice differentiable tp Q = f y l + 2 — k + J ib + Vi [ k + ibdz 

provided tb, tb z — > as z — > 00. Hence eq.(2.15) can V dx J J 

be written as oo 

~ + f (iak +t - k +xx + k +zz )ibdz, (2.23) 

J (iaF t -F xx + F zz )i;dz = 0, (2.17) { 

since = k +x + k +z , the total derivative in x. 

and therefore cq.(2.11) (for i = 1) is an operator iden- Hence if this equation is valid for all continuous ip, 

tity only if F (x, z; t) satisfies then 

iaFt — F xx + F zz = 0. (2.18) Vt (x; t) = -2^k+ = -2 (k +w + , (2.24) 
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(so that V\ is of degree zero), and k + (x,z;t) satisfies and so we choose 
iak +t - k +xx + k +zz + Vik + = 0. (2.25) 




On the other hand, we choose the undressed oper- 
ator Aq 2 ^ as follows 



(2.26) 



(where a and <5 are arbitrary real values), and 
then substitute it into eq.(2.12) : A^ 2 ) (I + $ + ) = 
(I + $+) A*, 2) , we obtain 

a \ 

S Vx - V 2 ( V + / fc+^^z 
8 J 

a 
+ / fc+ ( <5 | ip z dz 
5 




{i>x ~ Kip) 



a 
I- | 6 | / fc +x V(fe, 
<5 



(2.27) 



where we have written A' 2 ) = Aq 2 -* + V<i. On inte- 
grating by parts in the first term, eq.(2.27) becomes 

oo 

= V 2 jk+ 



+ 



OO / 

A- 

X \ 





OO / 


a 


+ 


/ 







X \ 





a 











5 











5 










6 





j ^ 





5 




a 











<5 











s 



a \ / a 

V2 = | (5 \k+ - k+ (5 
6 \ 6 



where fc+ (x, z; t) must satisfy 

k+ x +k +z 



(2.29) 




a 

(5 | +T/ 2 fc+ = 0. 
<5 / 

(2.30) 

We now return to the main development and de- 
scribe a further important step in the ZS dressing 
method. This is to introduce two pairs of operators 
AW and Aq\ A typical choice is in the following 
forms 

A^=Iia^--M; A< 2 > = 10-?- + L, (2.31a) 
at ay 

and 

A^ =Iia^-M 0] A< 2) - J/?|- + L 0) (2.31b) 

where /3 is constant, and M, L, M , and L are dif- 
ferential operators in x only. Consistent with our 
notation, Mq and Lq are comprised of constant co- 
efficients only and so A^\ A^ commute. Further- 
more, both Ag^and A^ are to commute with the 
same operator $i? so that eq.(2.11) is valid. On the 
other hand, operators A^ are defined according to 
eq.(2.12), with the same $+. Based on this point, it 
is instructive to examine the operator 



P = A«A( 2 ) (J + $+) - A^A« (J + $ H 
which, upon the use of eq.(2.12) twice, gives 



(2.32) 



= A«(/ + 



$+) A{, 2) - A< 2 ) (J + $+)A, 



A (1) A (2) _ 
A (D A (2) 



(I + $+)a( 2) a 



(i) 



(1) 



(2.33) 



However, A^and Aq 2 -* are chosen so that they com- 
mute with one another: hence P = 0. Thus wc obtain 



p = Ia^a^I (J + $+) = 



(2.34) 
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and, since (/ + <f> + ) is invcrtible (an easier assumtion 
which can be checked in specific cases), we get 



0. 



(2.35) 



which means that A' 1 ) commutes with A' 2 ). If we 
now introduce the choice given in eq. (2.31a) and 
eq. (2.316), eq.(2.35) becomes 

= (n4-M)(ipl y+ L 

-( //? i +L ) (^|- M >' (2 - 36) 



which simplifies to 



iaL t + /3M y + [L, M] = 0. 



(2.37) 



This is a generalisation of Lax pair to two auxiliarly 
variables; the Lax pair is recovered if we put (3 — 0. 
Equation (2.37) is the general system of the cou- 3.1 Type I 
pled NLS evolution equations which can be solved 
by the ZS dressing method (ZS scheme). We can 
also rewrite the connections among the operators L, 
M, M , L , AW and A( 2 ) related to the general cou- 
pled NLS equations as follows 



eq.(2.10)). This equations require F, and F is sup- 
plied by the solution of a pair of equations (eq.(2.18)), 
a pair since both equations (2.11) are to be satisfied. 
Note that the interesting point here is that the eigen- 
value does not appear explisitly at any stage in the 
ZS scheme. 

3 The New Integrable Coupled 
NLS Equations 

In this section, we investigate and derive two types of 
the integrable coupled NLS equations with their cer- 
tain perturbative terms and also show the Lax pairs 
of those equations using the ZS dressing method. For 
this purpose, we will initially choose certain matrix 
function k+, and then apply the ZS dressing method 
provided in section 2 to derive the equations we want. 



We start by choosing the following matrix function 
k + related to ZS dressing method as follows 



where 



and 



where 



A (2) 



L 



Vo. 



,(2) 



M = M -V 1 , 



d 2 



(2.38a) 



(2.38b) 



(2.38c) 



(2.38d) 



hi(x,z;t) i i 

h 2 (x,z;t) h s (x,z;t) 
±tF hi(x,z;t) h 5 (x,z;t) 



(3.1) 



hi, h 2 , h 3 , hi, and h 5 are functions which will be 
calculated using eq.(2.12). 

By substituting eq.(3.1) into eq.(2.24), we obtain 



Vi( X ,t) 



hi* 

K 



h2 x 
h 4 



w' 2 

h 5* J 



(3.2) 



On the other hand, by substituting eq.(3.1) into 
eq.(2.29), we find 



The procedure for solving eq.(2.37) can be de- 
scribed in the following explanations. The variable 
coefficients which arise in the 'dressed' operators, L 
and M, constitute the fuctions which satisfy the sys- 
tem of the evolution equations. These functions are 
known in terms of k + , where k + is a solution of the 
linear integral Marchenko equations (see eq.(2.9) and 



V 2 (x, t) = 



± 





u* 



(<*-&) 
u 






w 






(3.3) 



The functions hi x , h 2x , hs x , and in eq.(3.2) 
can be found by substituting eq.(3.3) and eq.(3.1) 



6 



into eq.(2.30), 

h\ x = T- 

h 3 x 



M 2 + H 2 ' 

\u\ \w\ 



and 



T 



5 — a 
S~ 

8 — a 
6 

5 — a 
S 

5 — a 



1 



u w 
1 

uw* 



w 



(3.4a) 

(3.4b) 
(3.4c) 
(3.4d) 

(3.4e) 



Hence, the coupled NLS equation type I (cq.(l.l)) 
and its complex conjugate can easily be derived by 
putting equations (2.21), (2.14), (2.26), (3.2), and 
(3.3) into eq.(2.35). 

Finally, we can obtain the Lax pair operators of 
cq.(l.l) and its complex conjugate using the relations 
of equations (2.37) and (2.38), 



a 





j dx 








S 
5 




/ 





(c-S) 


(a-S) 




± («-a) 
u* 


U 




W 




V 


± (S-a) 









,(3.5a) 



and 



M 



( h u 



dx 2 



+ 



h2 x 
h ± x 



h 3, 
/l5 



(3.5b) 



3.2 Type II 



Following the same procedure as type I, we only 
choose k + in a different form from that in eq.(3.1), 

ci(x, z,t) ^ 
k + = ( ^f- d(x,z;t) e(x,z;t) 

f(x,z;t) g(x,z;t) 

(3.6) 



a, d, e, / and g are are functions which 
will be investigated using eq.(2.12) and 

H (x, z; t) = ^1 + s(m 2 + H 2 ). 

By substituting eq.(3.6) into eq.(2.24), we obtain 



/ a 2 



Vt (x, t) 



( n ) fx 9- 



\ 



) 



(3.7) 

On the other hand, by substituting eq.(3.6) into 
cq.(2.29), we find 







V 2 {x,t) 



±(S-a)u q q 



±(6— a)w 







(3.8) 

The functions a x , d x , e x , fx and g x in eq.(3.7) can 
be found by substituting eq.(3.8) and eq.(3.6) into 
eq.(2.30), 



\u\ 2 + \w\ 2 



a 



1 + a (\u\ 2 + \w\ 2 Y 



d= + 



%= + 



S 

5 — a 
S 

S — a 



1 + 8 



(h 2 +h 2 ) 



U W 



l + s( 



u\ 2 + \w\ 2 



uw 



l + s(\u\ 2 + \w\ 2 ] 



and 



1 + 3 



(h 2 + h 2 ). 



(3.9a) 

(3.9b) 
(3.9c) 
(3.9d) 

(3.9e) 



Hence, the coupled NLS equation type II (cq.(1.3)) 
and its complex conjugate can be derived by sub- 
stituting equations (2.21), (2.14), (2.26), (3.7), and 
(3.8) into cq.(2.35). 
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We can then obtain the Lax pair operators of 
eq.(l.l) and its complex conjugate using the relations 
of equations (2.37) and (2.38), 




(3.10a) 



and 



M = I 



dx 2 



/ „V 



d x 

fx 



(f ), \ 



9x 



J 

(3.10b) 



4 The solutions of the inte- 
grable coupled NLS equa- 
tions 

In this section, we investigate the solutions of both 
coupled NLS equations (type I and type II). We 
firstly choose the general matrix function F, 

A(x,z;t) B(x,z;t) 

F = | A* (x, z; t) 

B*(x,z;t) 

'(4.1) 

where the functions A(x,z;t), B(x,z;t) and their 
complex conjugate will be found by solving the lin- 
ear differential equation in eq.(2.18). Hence, after 
substituting eq.(4.1) into eq.(2.18), we get the fol- 
lowing four linear differential equations 



iaA-t + A zz A xx 
iaB t +B ZZ - B xx 
iaA* t +A* z: 



iaBi + B* 



A* 
B* 



0, 

o, 
o, 

0. 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 



A 


= A e 


-apz 


e p(5x+ip(a 2 - 


s') t y 


B 


= B a e~ 


-apz 


e p(Sx+ip(a 2 - 




A* 


= A*e- 


-Scrz 


^a(^ax-\-ia(^5 2 - 


-a=")t)- 


B* 


= B*e 


-5<rz 


^cr^ax+ia^S 2 





The solution of the above equations can be derived 
by using separable variable method. We then find 

, (4.3a) 
, (4.3b) 
(4.3c) 
(4.3d) 

where Aq , Aq , Bq and Bq are arbitrary complex pa- 
rameters, and p and a are chosen as imaginary con- 
stants in which p* = —a. 

4.1 The solution of the type I 

Here, we derive the solution of the coupled NLS 
equation type I using the ZS dressing method ex- 
plained in section 2. By putting eq.(3.1) and eq.(4.1) 
into eq.(2.9), we get (for hi, u t and w) 

oo 

hl = _ f ^k e ^ 2 -« 2 Y e ~^ e <*°y dy 
J u 

X 

oo 

_y ^ e ^(^ 2 -« 2 )* e -^e^dj/, (4.4a) 



1 

u 



OO 

{Aq) J h 1 e i &( a '- s ') t e- a ""e Spu dy, 



(4.4b) 



and 



1 

w 



OO 

(B ) J h 1 e i &( a '- s ') t e- a ""e 



(4.4c) 
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The final solution is work on x = z as the ZS method 
works (see the explanation of eq.(2.11) in section 
2.2). Hence, by substituting eq.(4.4a) to eq.(4.46) and 
eq.(4.4c) we find the solution of the type I: 



u 



- (e- apz Ao) e 5px e l ^( a2 - s2 ) 

OO 



and 



where 



u (x, t) = — 
w (x, t) = — 



1 + 3DiD 2 

1 + SDiD 2 
B a 1 : 



and 



f(l^o| 2 + |B | 2 



((^)" 2 (P-)) 

©ii = e p( ~ 5 - a)x , 
e 2 i = e- CT(5 - Q)a: , 

e 12 = e-^(^- 2 >, 



e 22 = e ^- 2 )* 



□ l =011012, 



^2 — ©21 022- 



(4.5a) 
(4.5b) 

(4.6a) 

(4.6b) 
(4.6c) 

(4.6d) 

(4.6e) 
(AM) 

(4.6g) 



(4.7b) 



and 



= -K(e- apz Bo)e Spx e l£ ^( a2 - 52 ) t 

OO 

-N(B ) y ae i £(«'-*') t e- c ""e Spv dy. 



(4.7c) 

Hence, by substituting eq.(4.7a) to eq. (4. 7b) and 
eq.(4.7c) we find the solution of the type II: 



u (x, t) 
and 



-sTj 



l + ^□ 1 n 2 



(4.8a) 



«(..«> -I^Vr^^ I. <«»> 



where 



i -sT 2 y V 1 + ^ D i a 2 

|Ao| 2 DiDi 



Ti 



So, we can clearly see that the solution of the type I 
has no soliton solution. 



4.2 The Solution of the Type II 

By reviewing the solution's steps in section 4.1, we 
can also investigate the solution of the type II. We 
then start by substituting eq.(3.6) and eq.(4.1) into 
eq.(2.9), we get (for a, u. and w) 



oc 

J ^ L ^{S--^)t e -5a Zea av dy 

X 

OO 

_y ^ e ^ 2 -« 2 > e -^e^y, (4.7a) 



|B | 2 DiDi 



5 = 1 + 9fDin 2 + 9fnin 2 + 3 2 rJiD 2 nin 2 



Oi 



= 



Ti 



(4.9a) 
(4.9b) 

2> 

(4.9c) 
(4.9d) 

(4.9e) 



(l-I^)(l--T a )' 



(i-r^r)a-^i)' 

0' n = eP'(*-«)* } (4.Qf) 

21 = e -* (4.9g) 

6' 12 =e 4 ^( 52 -" 2 ) t , (4.9h) 
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e ^ = e-^(* 2 -« 2 >, (4.9i) 

°i = ©ii©^, (4.9j) 

and 

= e' 21 e 22 . (4.9k) 

It seems that the solution of the type II in equa- 
tions (4.8a) and (4.86) can roughly be identified as a 
stable soliton solution when parameter s is small. 
If we put s — ► either in eq.(1.6) or in equations 
(4.8) then we will find the one soliton solution of the 
coupled NLS equation of Manakov type appeared in 
ref.[6]. 

We can also choose parameters a = 4.5 and S = 
1.73 in our type II equation and its solutions, respec- 
tively. The aim of this choice is to make parameters 
X = \i = ^and to reduce our one soliton (in eq.(4.8a) 
and eq.(4.8&)) to be a simple solution. 

It is very interesting to expand the discussions re- 
lated to an equation which has a soliton solution. 
So, it is undoubtedly that the type II equation is 
particularly to be discussed in this paper. Even the 
equation (1.3) with a simple additional term has ex- 
perimentally been investigated by Christodoulides et 
al. [7], and Kivshar et al. [8]. Based on their ex- 
panations, we suppose that there must be a physical 
meaning of our type II equation closely related to 
that in ref.[7, 8] . 

5 Concluding Remarks 

In this paper, we have presented two types of the 
integrable coupled NLS equation with perturbative 
terms derived using the ZS dressing method and have 
also found their solutions. Using the Lax pairs of the 
two types, we have identified the existence of the in- 
tegrability of those equations. The very important 
step to obtain the two types of the coupled NLS 
equations in the ZS dressing method is the choice of 
k + in eq.(3.1) and eq.(3.6), respectively. The type I 
has no soliton solution. However, the type II could 
exactly have a stable soliton solution when the effec- 
tive saturation parameter, s is small (s — > 0). 

In the two types, the coupled NLS equations 
maybe have exclusively physical meanings. They 



may describe wave propagations in birefringent op- 
tical fibers with nonlinear effects such as the Kerr 
effect, The Raman Scattering and gain and loss ef- 
fects. 

The type II is very interesting to be discussed be- 
cause it has a complicated soliton solution . We also 
propose that this type must also have multisolitons 
solutions. However, the cases will be disscussed in 
our next investigations [10]. 
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